import sympy as sp
import math
import matplotlib.pyplot as plt
import numpy as np

# 画出速度曲线图
t = np.arange(1, 3400, 0.01)
v = 0.008

w = (1.08 + np.sqrt(1.1664 + 0.000064 * t)) / (2 * r)
xita = 1.08-w*t
print(xita)
vx = 50 * xita * np.cos(xita + 1.082) * (-w)
plt.figure(figsize=(6, 6))
plt.plot(t, w)
plt.plot(t, vx)
# plt.ylim(-2.5, 30)
# plt.xlim(0.1, 3500)
plt.show()
